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EXPLORING LIMIT BEHAVIOUR OF NON-QUADRATIC 
TERMS VIA H-MEASURES. APPLICATION TO SMALL 
AMPLITUDE HOMOGENISATION. 

MARTIN LAZAR 


Abstract. A method is developed for analysing asymptotic behaviour of 
terms involving an arbitrary integer order powers of functions by means 
of H-measures. It is applied to the small amplitude homogenisation problem 
for a stationary diffusion equation, in which coefficients are assumed to be 
analytic perturbations of a constant, enabling formulae for higher order correc¬ 
tion terms in a general, non-periodic setting. Explicit expressions in terms of 
Fourier coefficients are obtained under periodicity assumption. The method 
enables its generalisation and application to the corresponding non-stationary 
equation, as well as to some other small amplitude homogenisation problems. 


1. Introduction and motivation 

H-measures, as originally introduced a quarter of century ago by L. Tartar [17] 
and (independently) P. Gerard [10] are kind of a microlocal defect tool, measur¬ 
ing deflection of weak from strong L^ convergence. They explore quadratic limit 
behaviour of bounded L^ sequences. 

A prominent feature of H-measures is their capability to keep track of an equation 
satisfied by functions generating them. More precisely, if an H-measure is associated 
to solutions of an equation = 0 (accompanied by a series of initial/boundary 
conditions), one can take advantage of their basic properties: the localisation and 
the propagation one. The former constrains the support of the H-measure within 
the characteristic set of the (pseudo) differential operator P, while the latter states 
that the measure (as well as concentration and oscillation effects) propagates along 
bicharacteristics of P. 

Since their introduction, they have been successfully applied in many mathemat¬ 
ical fields - let us here mention generalisation of compensated compactness results 
to equations with variable coefficients [10, 17], applications in the control theory 
[9, 16], the velocity averaging results [10, 13], as well as explicit formul® and bounds 
in homogenisation [17, 4, 5]. 

Original H-measures are restricted to the L^ framework. This constraint has 
partially been overcome with the introduction of H-distributions [7] - a generali¬ 
sation of the concept to the L^’,p > 1 framework. More precisely, a new tool is 
constructed with the aim of exploring products of a form f UnVn, with a sequence 
(m„) being bounded in L^, while Vn are taken from the corresponding dual IP 
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However, the above generalisation is still bound to a study of quadratic terms 
(with possibly different factors), and does not handle higher order terms, such 
as cubic ones. The aim of this paper is to investigate possibilities of handling a 
general LP,p > 2 sequences and to describe, roughly speaking, a limit behaviour 
of / \un\^- Here it should be mentioned that actual products we deal with are 
more complicated, allowing each factor to be accompanied by a pseudodifferential 
operator acting on it. 

The analysis below is based on a relatively simple idea: for a sequence of functions 
Un converging weakly to zero in L^’,p > 2, the sequence is bounded in L^, and 

determines (up to a subsequence) an H-measure denoted by ^p/ 2 - Thus the limit 

of / \un\'^ = f can be expressed by the measure /ip/ 2 , which in some cases 

can be calculated explicitly. Of course, one has to be aware that in the absence 
of strong convergence, the weak limit of in general differs from zero, which 
requires correction terms entering the expression. 

The obtained results are applied to a small amplitude homogenisation problem. 
The main idea of a small amplitude approximation consists of taking a (formal) 
expansion of a solution to a problem under consideration with respect to a small 
parameter representing perturbation of the coefficient. Originally introduced by L. 
Tartar in [17] for a stationary diffusion problem, the approach has subsequently 
been elaborated and applied to more general homogenisation [4, 8] , optimal design 
[1, 2] and inverse problems [12]. In all these papers H-measures are used as the 
main analytical tool, and the analysis is performed up to the second order expan¬ 
sion, within which the result is obtained by exploring limits of quadratic terms. 
Handling of non-quadratic terms, appearing in higher order expansion, is, how¬ 
ever, not achievable by a standard usage of H-measures and requires a different 
approach. In this paper we try to make a step forward in that direction and to 
obtain expressions for higher order correction terms. 

The paper is organised as follows. In the next section we describe in details the 
above presented idea, and demonstrate the method of expressing limit of higher 
order terms via (original) H-measures associated to an appropriate combination of 
given LP functions. Application of the method is demonstrated in Section 3 on the 
example of small amplitude homogenisation for a stationary diffusion problem, with 
a particular intention given to a special, but important case of periodic functions. 
The paper is closed by some concluding remarks, and by pointing toward some 
related and open problems. 


2. Exploring non-quadratic terms through H-measures 

2.1. General setting. The original H-measures explore quadratic limit behaviour 
of terms, and their existence theorem can be formulated as follows (cf. [10, 17]). 


Theorem 1. Let {un) be a sequence converging weakly to zero in L^(R'^). Then, 
after passing to a subsequence (not relabeled), there exists a nonnegative Radon 
measure fi on the cospherical bundle x sueh that for every (pi, p 2 G Co(R'^) 
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and Ip G C(S“ it holds 


lim 

n—voD 


.4^(<PlUn)(x)((/32Un)(x)dx = (/i, ((^ 1 (^ 2 ) ^ 4 ’) 


( 1 ) 




(/3i(x)(p2(x)'0(^)d/i(x,^) . 


R‘ixS<i-i 

where A,p is the (Fourier) multiplier operator on associated to 

The measure p, is called the H-measure associated to the (sub)sequence (u„/). 
In the sequel we shall often abused terminology and notation by assuming that we 
have already passed to a subsequence determining an H-measure. 

Notation. Throughout the paper by (•, •) stands for a sesquilinear dual product, 
taken to be antilinear in the first, while linear in the second variable. By Kl we 
denote the tensor product of functions in different variables. □ 

The theorem also generalises to the vector sequences (un), resulting in a hermit- 
ian, positive semi-definite matrix H-measure, whose diagonal elements are (scalar) 
H-measures associated to a corresponding component ul^. 

The proof of the last theorem is based on the (First) commutation lemma [17, 
Lemma 1.7], enabling exchange of the multiplier operator and the operator of 
multiplication by Lpi in (1) when passing to the limit. Thus the limit depends on 
the product ipi(p 2 only and results in a bilinear functional on Co(R‘^) 0 C(S‘^“^). 
Taking into account its positivity, and by using the Schwartz kernel theorem, as well 
as the Schwartz lemma on nonnegative distributions, one shows that the functional 
is a Radon measure in both variables, x and 

A multiplier operator associated to a bounded symbol ip is a continuous op¬ 
erator on L^(R‘^), which is easily demonstrated by means of the Fourier transform. 
Generalisation of that result to an setting is provided by the Marcinkiewicz mul¬ 
tiplier theorem (e.g. [11, Theorem 5.2.4] showing that A^ is a bounded operator 
on LP(R‘^), for any p G (1, oo) and ip of class C'^. 

We would like to generalise Theorem 1 by considering higher order expressions 
in G LP,p € N. More precisely we are interested in expressing the limit of 


/R<i 


A-4,^{(piUn)(x.)A^^{ip2Un){:>i) ■.. A^^{(ppU„){x)dx, 


( 2 ) 


where (pi^ipi,i = l..p are appropriate test functions. The following result holds. 

Theorem 2. Let (u„) be a sequence converging weakly to zero in L^^'^(R'^) for 
some p G N and e > 0. Then for any choice of test functions (pi G Cc(R‘^),ipi G 
C^(S^-i),i = l..p it holds 

lim/ A- 4 ,^{ipiUn){x) ■ ... ■ Aiij„{ippUn){x)dx = {pyy,,(pMl) + {ipv){x)w(x)dx, 

^ dRd J 

( 3 ) 

where ip = Y\d=i 7’u while is the off-diagonal component of the matrix H- 
measure associated to the sequence (u„ — v,Wn — w) defined by (5) and (6) below. 

Proof: By the theorem assumptions, the products (fiUn are bounded in LP+®(R'^) 
for all i = l..p. Additionally the Marcinkiewicz multiplier theorem provides bound¬ 
edness of AtPi {(piUn) in the same space, as well. As by the (commutation ) Lemma 3 
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the commutator Ci = Aip^ipi — is a compact operator on L^’(R'^), we can con¬ 

secutively exchange the order of operators in (3), showing that the corresponding 
limit equals 


lim / ipA^^{xiUn){-x.) ■ ... • A^(xpUn){x)dx = lim / ipVnWndx. 

n n J^d 


where the last two factors are defined by 


Vn ■ = 


(p-l)/2 


p/2 

Y\ AtiiXxiUn), p e 2N 

i=l 


/ xl/2 

n A^,(xiU„)IA^._^_^.,Jx(p+i)/ 2 U„) , elsewhere 

i=i ^ ^ 


and 


Wn : = 


fl Aji.iixi'Un), p G 2N 

i=p/2-|-l 

/ \l/2 P 

{■^■ 4 >(p+i)/ 2 iX{p+i)/ 2 Un)) n A^^ixiUn), elsewhere. 

^ ^ i={p+3)/2 


(4) 


(5) 


( 6 ) 


In the above Xi stands for the characteristic function of the support of ifi. The 
conjugation sign originate from the scalar product used in the definition of H- 
measures (1). We also assume that one branch of the square root has been selected 
(the choice does not effect the result). 

Sequences (u„) and (w„) are bounded in L^(R‘^), and (after possible passing to a 
subsequence) we denote their weak limits by v and w, respectively. Note that in 
the absence of the strong (zero) convergence of (u„), the limits v and w in general 
differ from zero. 

Rewriting the last term in (4) as 


lim 

n 



<f{Vn 


v)(Wn 


w)dx + 



ipvwd'K, 


and expressing the last limit via H-measure one obtains the result. □ 

The generalisation of the First commutation lemma to the non-L^ framework 
can be found in [7], whose (slightly corrected) version we reproduce here. 

Lemma 3. Let {un) be a bounded sequence of functions in L^(R^) n L^’(R‘^) for 
some p G (2,oo], converging weakly to zero (in the sense of distributions). Denote 
by C = Axpip — pA.^ the commutator determined by ip € Co(R ^) andfjG C'^(S'^-i). 
Then the sequence (Cun) converges to zero strongly in L'J(R'^) for any q G [2,p). 

The proof of the lemma is based on the interpolation inequality, the First com¬ 
mutation lemma ([17, Lemma 1.7], providing compactness of the commutator C on 
L2(R'i)), and the Marcinkiewicz multiplier theorem (providing boundedness of the 
multiplier A.^ on L'J(R‘^) for any g G (1, oo) and i) smooth enough). 

Remark 4. • The additional regularity assumptions on (m„) (requiring it to 

belong to L^’+®(R'^ ) for some positive e) and ip (assuming it to be of class 
C"^) are necessary for the application of the (commutation) Lemma 3. It 
enables the limit in (3) to depend on the product p of test functions pi 
only. 
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• Note that by relation (4) it follows that the limit in (3) defines a continuous 

(p+ l)-linear form B on Co(R‘^) x x ... x In the special 

case of p = 2 one can use Plancherel theorem in order to get a bilinear 
form depending on products p = ipi and only, resulting in 

a distribution, and eventually an (H-)measure on R'^ X Usually this 

tool is not at our disposal, and we are left with a more general object. 
However, relation (3) shows that the object is a measure with respect to x 
variable. More precisely, for any choice of test functions = l-.p, the 
functional equals the measure +vwX{x), with the last 

term standing for the Lebesgue measure. 

• The last theorem can be adapted in order to compute limits of more gen¬ 
eralised expressions of the form (for simplicity of the notation we restrict 
to the even p): 

j • . . . • (‘^p/2'Wn)^ (x) (7) 

• ■^V’p/2+i(V5p/2+iw») • • ■ • • Mv,,p(ppU„)(x)dx. 

The difference compared to (3) is in an additional multiplier (with symbol 
tjj) acting on half of the factors. 

Carefully repeating steps of the proof of Theorem 2 one obtains that the 
limit of (7) equals 

(Ppu, , p KI V') + J ...■ Pp /2 v) (x) ipp/ 2 +l ■ ■■■■ (Ppw{-x)d'x. 

the result coinciding with (3) for i/; = 1. 

2.2. Localisation principle. As already stated in the introduction, one of the 
most important properties of H-measures is the so called localisation property (e.g. 
[10, Corollary 2.2]), playing a crucial role in most of successful applications of the 
tool. Thus it would be important to check if the same principle applies when 
analysing limit behaviour of higher order terms. Here we substantiate a positive 
answer for an L°° sequence, but it holds for more general situations as well. 

Suppose Un & L°°(R‘^) are solutions to an equation = 0, converging vaguely 
to zero. By the localisation principle, the associated H-measure is supported within 
the characteristic set of the operator P. But the same constraint also applies for 
off-diagonal terms of a matrix H-measure associated to a vector sequence (u„, Wn), 
with (wn) being an arbitrary bounded L? sequence. 

For that reason, we can split the integrand in (2) into two parts, the first one 
coinciding with the first factor, while the second one containing the rest of them, 

p 

i.e. Wn = n A,p-{(piUn)- Taking into account Remark 4 (the last item), we get 

i=1 

the limit of (2) equals {puw,^ where puw is an off-diagonal term of an H- 

measure associated to (u„,w„ — w). In accordance to above, it satisfies the same 
constraints as an H-measure associated to (un), i.e. its support is localised within 
the characteristic set of P. 

2.3. Explicit formulae in the periodic setting. Let (m„) be a bounded sequence 
of periodic functions in Lj^j,(R‘^) defined by 

u„(x) = ^ 
keZ'* 


(8) 
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Hereby we assume the mean value uq equals zero, thus obtaining the zero weak 
convergence of (u„). The associated H-measure is well known, and is combination 
of the Lebesgue measure (in x) and the Dirac mass (in ^): 

k (l)A(x). (9) 

|k| 

k 

By using the above described procedure, we would like to express the limit of (2) 
in the case p = 4. First, let us note that for a periodic function the action of 
the multiplier results in a periodic function again: 

^^w„(x) = ^iik^(k)e“-’^, 

k 

where we have used that the symbol of the above multiplier is a homogeneous 
function of order zero. Specially, function defined by (5) in this case takes the 
form 


U„(x) = {Atl,^Un){A^^Un){x.) = ^ Uj Uk V'l (j)V '2 (k) 

jF 

and converges weakly to 

7 ;(x) ^ Wk^-k' 0 i(k)V^ 2 (-k) . 

k 

Similarly, 

Zl^(x) = (^V’3^n)(^V’4^n)(x) = ^ 7/lWm'03(l)V'4(m) 

l,m 

1 

The measure pvuj determined by the sequences (u„ — v) and {wn — w) then reads 
/i™ = X! ( X! MjUkUlMmV’l(j)'</'2(k)l/l3(l)V’4(m)^ 5.^(^)A(x). 

j.k l,m 

j+k#{0} l+m=-(j+k) 


Taking into account the form of the limits v and w, relation (3) implies 


lim 

n 



A.4,^{(piUn){x) ■ ■ A^^{(p4Un){:ic)dx 


Wj'UkUlWm'*/'l(j)'*/'2(k)V'3(l)V'4(m) 

j,k,l,m 

l+m=-(j+k) 



V?(x)dx, 


( 10 ) 


where, let it be repeated, (p = 

Similar approach can be applied for other values of p. Here we briefly comment 
the case p = 3, as a prototype for odd values of p. 

Note that due to additional regularity of sequence («„) (bounded in LJ^^(R‘^), 
and thus in L^^^(R‘^), for any q), we can avoid splitting of the middle factor in ( 2 ) 
into two terms involving the square root, and define 

Tn — Aqp^Un^ Wji — (^^244n) Un) 5 


EXPLORING LIMIT BEHAVIOUR OF NON-QUADRATIC TERMS VIA H-MEASURES 7 


which both belong to Ljq^(R'^). Thus the associated measure is well defined and 
computing as above one gets 

fj-vw = X! ( X! WkWlMm'0l(k)'!/'2(l)t/'3(m)^ (5^(1) A(x). (11) 

k 5 ^{ 0 } l.m 

k+l+m=0 


and 


lim 

n 



A^^{ipiUn){y^) ■ ■■■■ A^s{ip3Un){x)dx 


= E 

k,l.m 


Uk ^1 -!/'i(k)V’2(l)t/’3(m) 



>p{x)dx, 


k+l+m=0 


( 12 ) 


Similarly, the last formula is now easily extend to an arbitrary integer p > 2, 
providing the following result. 

Theorem 5. Let (un) be a bounded sequence of periodic functions in 
defined by (8) , with the mean value zero. Then for any p G N and any choice of 
test functions (pi € Cc(R‘^), V”! G C'^(S‘^“^),7 = l..p it holds 


lim 

n 



Alv,i(piWn)(x)-. . .■A^^{(ppUn){x)dx = 


E (n^k,?/'i(k*) 


kiSZ" 


Eiki=0 


where ip = OLi Ti- 


R'i 


p(x)dx, 


(13) 


Remark 6. The last theorem easily generalises to a case when each factor of the in¬ 
tegrand in (13) is associated to a different sequence {u\f)n, i = l-.p, just by adjusting 
the Fourier coefficients on the right hand side. 

Note that the last theorem also incorporates the expression (9) for an H-measure 
associated to a sequence of periodic functions. 


3. Computation of higher order correction terms in small amplitude 

HOMOGENISATION 

3.1. Small amplitude homogenisation. We consider a sequence of elliptic prob¬ 
lems: 

J -div (A"Vm”) = / G H-^(fl) 

\ u" G H(,(r!), 

where 11 C R'^ is an open, bounded domain. 

If the coefficients A” are taken from the set (with 0 < a < (3) 

M{a,f3;^) ■= {A G L°“(ll; £(R^; R^)) : A(x)|.| > & A-i(x)|.| > (a.e. x G H)} , 

we have that (after extracting a non-relabelled subsequence) A"—^A°° G At (a, /3; LI) 

(e.g. [18, Theorem 6.5]). The above H-convergence means that for any choice of 
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/ G H the sequence of solutions to (14) satisfies: 

m” —^ in Hj(rj) 

^ in L 2 ( 11 ) , 

where is the solution of (14) with oo instead of n. 

For small amplitude homogenisation we consider A" to be a perturbation of a 
constant: 

A”(t,x) = Ao + , (15) 

where A" —^ 0 in L°“(n) for any i >1. 

Assuming that Aq G Ai{a, /3; fl), we have that (for small values of 7 ) 

A;-^A^ = Ao + 7 A^(t, x) + 7 "A“(<, x) + 7 ^A^(t, x) + 0 ( 7 ^) , (16) 

where the homogenised limit A“ is measurable in x and analytic in 7 (for details 
on small amplitude homogenisation consult [18, Chapter 29]). 

The components of the limit A“ up to the second order of 7 are given explicitly 
by means of H-measures. More precisely, the following theorem holds (cf. [17, 
Theorem 4.2]). 

Theorem 7. The effective diffusion tensor Aff satisfies 
A-=Ao + 7"A“(t,x) + o(7"), 
where the second order correction is given by 

[ (A^),,w<(>(x)dx = 

with fill standing for an H-measure (with four indices) associated to A". 

3.2. Expressing higher order correction terms. The goal of the following 
study is to use the approach presented in Section 2 in order to construct explicit 
expressions for higher order homogenised terms. We shall see that such obtained 
formulas will include H-measures associated to (the powers of) the lower order per¬ 
turbations. Here we follow the approach presented firstly in [17] for elliptic problems 
(for the parabolic version see [4]). 

Fix a function u G HQ(r2), and denote by m” the solution of 

-div(A”VM”) = -div(A“Vu). (17) 

Because of H-convergence, we have that ut) —^ u (in H(j(r2)) and that D" := 
A”Vm" —^ AffVu (in L^(H)). After writing the expansions in powers of 7 : 

u; = u^+ 7 < + 7 ^ul) + 0 ( 7 ") , D!) = + 7 D- + 7 ^ 0 ^ + 0 ( 7 ^) , 

we see that Uq —^ u and u" —^ 0 for i > 1 , while it remains to calculate the limits 
of 

Equating the terms with equal powers of 7 we get Vuq = Vm, Dq = AqVm (with 
7 °); and D" = AqVu" -I-A”Vm —^ 0 (with 7 ^). On the other hand, D" converges 
to A“Vu (the term in expansion of AffVu with power of 7 equal 1 ), which, by 
varying u, gives that A“ = 0 . 

For the quadratic term we have: 

D” = AqVm^ -f A”V< -f A”Vm ^ lim A”V< = A“Vm , 

n 


(18) 
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By expansion of (17) and taking terms with 7 ^, observe that u" satisfies the problem 
(14) with A” = Aq and the right hand side equal to div(A"VM). As Aq is a 
constant matrix, application of the Fourier transform yields 


Vu”(4) = 


(|®g) (A"Vu)(|) 

Aol-I 


or equivalently 

V<(x) = -A^ (A^Vu) (x), (19) 

where is the multiplier with the symbol 41(4) = 

Taking into account (18), one gets that the second order correction A 2 is ex¬ 
pressed via limit of the quadratic term in A", i.e. via a corresponding H-measure, 
which is essentially the result of Theorem 7. 

In order to find the cubic term, note that 


D” = AoVu” + A”Vu” + A^V< + A^Vu lim(A”VM” + A^Vu”) = A“Vu , 

( 20 ) 

where we have taken into account zero weak convergence of (A") and (u") for i > 1. 

On the other hand, by equating the terms with 7 ^ in the expansion of the relation 
(17), one gets for U 2 

-div (AoVu^) = div (A^Vm + A^Vu^ - A“Vu) ^ 0, 

where the convergence of the right hand side follows from (18). 

By applying the Fourier transform, similarly as in (19) one gets 

Vu^(x) = -Avp (A^Vu + A^V< - A“Vu) (x). 

Putting the last expression together with (19) in (20), integration and multiplication 
by a test function (p G C(f(fl) yields 

J (^A“Vudx = lim J ip(^-A^A^A^Vu-Al^A^A^S/u+A^A^ (A^A^pA^Vu))dx 

( 21 ) 

The limits of the first two terms on the right hand side are expressed via H-measures 
determined by functions A” and A^- The last term in (21) is a cubic term in A". As 
applications of H-measures so far have been constrained to quadratic expressions, 
for that reason no expression for the third order correction has been obtained. In 
the sequel we shall overcome the restriction and calculate its limit by means of 
Theorem 2. 

We analyse the quadratic terms first. Carefully expanding the product by com¬ 
ponents we obtain 

lim / pfA^A^A^ + A^A^A^^Wudx = (25R/Xi2, , 

^ Jn ^ ^ Ao4-4 

where fii 2 is an off-diagonal block term of the H-measure associated to (A", A 2 ) - 
a measure with four indices (the first of them not being contracted above). 

As for the cubic term, observe that by using an pseudodifferential calculus iden¬ 
tity f (At/ju) vdx = j uA^vdx (where ip denotes the change of the argument sign, 
i.e. i/>( 4 ) = ip{—0) it can be rewritten as 

\imj (^A^iipA^fY A^A^A^VudTC, 


(22) 
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where symmetry and evenness of has been taken into account. Denoting 

V" = ^ , W” = (23) 

Theorem 2 gives that the i-th component of (22) equals ^ ‘fdju 0 1) 

denotes appropriate component of the four-index H-measure /Xyn/ asso¬ 
ciated to (V",W )). Note that we omit the last term appearing in (3), as the 
(zero) weak convergence assumption on coefficients A" implies the weak limit of 
V” equals zero. By varying function u S C^(D) (e.g. choosing Vu constant on w, 
where a; (e D) we finally obtain the following result. 

Theorem 8. The third order correction of the effective diffusion tensor defined 
in (16) is given by 

[ i-^Tr Tdyi = , (24) 

Jq, -^-04 • 4 

where denotes the matrix H-measure determined by sequences (A") and (A^) 
with components , and similarly for the H-measure ffyyy, whose 

generating sequences are given by (23). 

Remark 9. Although not given explicitly, the last term in (24) depends on ik = 
indeed, as both sequences (D"),(1T”), as well as the associated H-measure 
depend on it. This dependence is more apparent if a periodic setting is considered, 
which is the subject of the next subsection. □ 

Similar procedure as in the above can be applied in order to express further order 
correction terms of the effective diffusion tensor. Derivation for the fourth order 
term is sketched in the next subsection under periodicity assumption. 

3.3. Periodic setting. This subsection is devoted to a special, but important 
case of periodic coefficients in (14). By employing techniques of Section 2 we derive 
explicit formul® for higher order terms expressed by Fourier coefficients of A". 

We suppose the coefficients in (15) to be of the form 

A”(nx) = A,(nx) = ^ i € N 

IcGZ-i 

where (just to simplify calculations) each A^ is symmetric with the zero mean value. 
The measure associated to (A", A 2 ) takes the form 

/X 12 = ^ Akk ® A 2 _k(5^ (l)A(x). 

k 

Thus the first term on the right hand side of (24) equals 

-2y^ I / (Ai,kk) 0 (A 2 _kk) (^(x)dx. 

^ Aok ■ k Jq 

As for the second term one uses the formula (11) for an H-measure associated to 
cubic terms (in A" here). Thus the trace of (^ke trace being taken with 
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respect to the middle two components) equals 

A .1, ™ 



k.l.m 

k+l+m=0 


and for the third order correction we obtain the explicit expression 



k+l+m=0 


Just for comparison, here we provide an analogous expression for obtained by 
means of Theorem 7 : 



By virtue of Theorem 5 derivation of expressions for higher order terms goes sim¬ 
ilarly, just involving more lengthy computations of linear algebra. Here we briefly 
provide the result for A^. 

Similarly as (21) one obtains 




A^Vudx = lim / ^ A^A^,A^Vu - A^A^,A^Vu - A^A^A^Vu 

n In V 


+ (A^,A^)^ A'^A^A^Vu + (A>kA”)^ A^A^A^^Vu 

+ (A^,A”)^ A^A^A^Vu - (Am,W")^ W”Vu)dx 


where W" is given by (23). 

The first three terms on the right hand side are quadratic terms whose limit is 
expressed via standard H-measures. The next three are cubic ones, and the limit 
is obtained by virtue of relation (12), similarly as it was done above for A^. The 
last one is a quartic term in A” (as W" are quadratic expressions of A"), which is 
treated by means of formula (10). Performing some tedious, but mostly elementary 
computations one gets 



k+l+m=0 


+ Ai,ik • m (A 2 ^kk) (g) (Ai,mm) 


l+m=-(j-(-k) 
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Remark 10. • The above results are easily generalised to the case when 

Af - ^ Ai ^ 0, with correction terms entering expression for the ho¬ 

mogenised limit A“. 

• In analysis of small amplitude homogenisation we considered diffusion co¬ 
efficients A„ to be perturbations of a constant. This was important in 
order to explicitly express Vu" via a given function u (relation (19) etc). 
However, one can try to generalise result to a variable Ag by justifying and 
following suggestions given in [8] or [17], partially based on the localisation 
principle for H-measures. Here it is important that the measures analysed 
in previous section obey the same principle in accordance to subsection 2.2, 
making potential generalisation feasible. 

4. Conclusion 

In the paper we have developed a method for expressing limits of non-quadratic 
terms by means of original H-measures, followed by application to the small am¬ 
plitude homogenisation problem for a stationary diffusion equation. The method 
provides higher order correction terms of the effective diffusion tensor expressed by 
H-measures associated to non-quadratic terms in a general, non-periodic setting. 
Explicit formul® in terms of Fourier coefhcients are obtained under periodicity 
assumption. 

A similar approach can be also conducted in the case of non-stationary diffusion 
problems by means of parabolic H-measures [4, 6] . The latter represent a generali¬ 
sation of the tool that takes into account the difference between the time and space 
variables, intrinsic to parabolic type problems. 

The corresponding analysis, as well as the explicit formula for the second order 
homogenised term has been provided firstly by [4] for the periodic functions, while 
non-periodic generalisations are given in [8] . Results and the approach presented in 
Section 2 of expressing non-quadratic limits via microlocal defect measures extend 
easily to parabolic H-measures as well. Thus the third order homogenised terms 
for a non-stationary diffusion problem can be easily obtained by slightly adapting 
the proof of Theorem 8, and similarly for higher orders. 

The presented method for expressing non-quadratic terms can also be generalised 
to recently introduced multiscale H-measures [19, 3]. Unlike the original ones, mul- 
tiscale variant are capable of distinguishing sequences with different frequencies. 
Thus, their application in expressing higher order terms in small amplitude ho¬ 
mogenisation makes full sense in the case of coefficients A" oscillating on various 
scales. 

The method developed in this paper paves the path for more precise analysis of 
other problems involving small amplitude homogenisation [1, 2, 12] in which the 
analysis so far has been conducted by means of H-measures up to the second order 
expansion. It is important to notice that in all these papers sequences describing the 
limit terms posses higher regularity than the one required by the original definition 
of H-measures (L°“ instead of merely one) and fit within the setting of the 
presented approach, thus making its potential applications feasible. 
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